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Aiming at the optimization of the spin diffusion length in (001) GaAs quantum wells, we explore 
the effect of the anisotropy of the spin-orbit coupling on the competition between the Rashba and 
the Dresselhaus spin-orbit couplings by solving the kinetic spin Bloch equations with the electron- 
phonon and the electron-electron scattering explicitly included. For identical strengths of the Rashba 
and the Dresselhaus spin-orbit couplings, the spin diffusion length shows strong anisotropy not only 
for the spin polarization direction but also for the spin diffusion direction. Two special directions are 
used seeking for the large diffusion length: (110) and (110). Without the cubic term of the Dressel- 
haus spin-orbit coupling and with the identical Dresselhaus and Rashba strengths, infinite diffusion 
lengths can be obtained either for the spin diffusion/injection direction along (110), regardless of 
the direction of spin polarization, or for the spin polarization along (110), regardless of the direction 
of the spin diffusion/injection. However, the cubic Dresselhaus term cannot be neglected, resulting 
in a finite spin diffusion length which decreases with the temperature and the electron density. The 
anisotropy for the spin diffusion direction and spin polarization direction is maintained. For the spin 
diffusion/injection direction along (110), the spin diffusion length increases first with the increase 
of the Rashba strength (from 0) which can be tuned by the external gate voltage; when the Rashba 
strength is slightly smaller than (in stead of equal to) the Dresselhaus strength, the diffusion length 
reaches its maximum, followed by a decrease with further increase of the Rashba strength. 

PACS numbers: 72.25.Rb, 72.25.Dc, 71.70.Ej, 71.10.-w 



I. INTRODUCTION 



The control of the electron spin in nanoscale is an ob- 
ject of great interest in the field of spintronicsJ^ The 
knowledge of how the spin evolves as the electron goes 
through a nano-device is a key element for the purpose 
of using the spin degree of freedom as a mechanism for 
information transfer and processing. Randomization of 
the electron spin, however, occurs due to the spin-orbit 
coupling as the electron wave vector changes during and 
after scattering processes. Between scattering events, the 
precession frequency of the spin about the local magnetic 
field also changes from electron to electron depending on 
the wave vector. To overcome these detrimental effects, 
many efforts have been made to understand the spin dif- 
fusion and the spin relaxation^&l&^ill^i^i^ 
in particular after the proposal of the spin-field-effect- 
transistor (SFET) by Datta and Das more than sixteen 
years agOfil In such a device, ferromagnetic material is 
used as source and drain for the injection and the de- 
tection of spin-polarized electrons. The electron spins, 
moving through a quasi-one-dimensional channel, pre- 
cess about an effective magnetic field tuned by a gate 
voltage. In n-type zinc-blende semiconductors the lo- 
cal magnetic field gives origin to a Zeeman-like splitting 
which, combined with the scattering, can cause a spin re- 
laxation/dephasing known as the D'yakonov-Perel' (DP) 
mechanism.— In GaAs, the Dresselhaus terrn^ is dom- 
inant in the DP wave vector dependent magnetic field, 
giving rise to a spin splitting due to bulk inversion asym- 



metry. Furthermore in low-dimensional semiconductor 
structures with asymmetric confining potential, another 
term known as the Rashba term2£ contributes to the DP 
mechanism, giving rise to a spin splitting based on the 
structure inversion asymmetry. While the Rashba term 
is linear, the Dresselhaus term is cubic in the wave vector 
components. 

A few years ago Schliemann et al^ proposed a non- 
ballistic spin-field-effect transistor based on the competi- 
tion of the Dresselhaus and the Rashba terms. In such a 
transistor a gate voltage is tuned to give equal strengths 
to both terms, leading to a very long spin dephasing time 
for the spin polarization along (110) direction^ In their 
work the cubic term is argued to be unimportant. Cheng 
and Wu studied the effect of the cubic term on the spin 
relaxation time of spins along the (110) direction by solv- 
ing the kinetic spin Bloch equations, and obtained a fi- 
nite spin relaxation timei^l Actually, in spin transport, 
when spins are injected into a sample, in addition to the 
direction of spin polarization, the direction of spin dif- 
fusion/injection is also important. The spin relaxation 
in transport depends on both directions. In the present 
work, the anisotropy of the spin-orbit coupling is taken 
into account when we study the ideal balance between 
the Rashba and the Dresselhaus terms in a quasi-two- 
dimensional channel of a GaAs quantum well. Therefore, 
both the direction of the spin polarization and the direc- 
tion of the spin injection are analyzed. We further show 
that the cubic term appearing in the Dresselhaus effective 
magnetic field cannot be neglected in the competition 
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FIG. 1: (Color online) Schematic of the different directions 
considered for the spin polarizations [(110), (110) and (001)- 
axes] and spin diffusion/injection (z-axis). 



with the Rashba field. Our treatment goes beyond the 
single particle method in that we use a fully microscopic 
many-body approach by solving the kinetic spin Bloch 
equations (KSBE). 4 As shown by Weng and Wu, the cor- 
relations between the spin- up and spin-down states, i.e., 
the off-diagonal terms of the density matrix in spin space, 
play an essential role in the spin diffusion/transport. 5 
Studying the spin polarizations and the oscillations in 
the direction of diffusion one can calculate the injection 
length as a function of the direction of injection. 

This paper is organized as follows: In Sec. Hi] we de- 
scribe the Rashba and the Dresselhaus spin-orbit cou- 
plings and study the effective magnetic field under differ- 
ent injection directions in the quantum well. In Sec. IIIII 
we construct the KSBE and we discuss how they change 
with the injection direction. In Sec. IIVI we apply the 
KSBE to our model and present the results for different 
parameters. We conclude in Sec. [V] 



II. SPIN-ORBIT COUPLING 

The DP mechanism introduces a spin-orbit term in 
the Hamiltonian which can be expressed by a Zeeman- 
like term with a wave-vector-dependent effective mag- 
netic field: 



For convenience, the effective g-factor and the Bohr mag- 
neton are absorbed into the definition of f2(k). For 
GaAs, which lacks an inversion symmetry center, the ef- 
fective magnetic field contains the Dresselhaus term. For 
electrons in infinite-barrier-height quantum wells with a 
small well width a in the Tq band, this term is;-- 



M*3 - (k 2 z )) 

fi„(k)= 7 | ky{{k*)-& 





(2) 



where 7 represents the spin splitting parameter and 
{kl) = {v/aY. 

If a space gradient is applied to the sample, another 
effect shows up, the Rashba effect. The Rashba term in 
the Hamiltonian for electrons in the T§ band is 



H R = 



r\f c a • k x e, 



(3) 



where r|j 6c is a material parameter, and e is the electric 
field determined by the asymmetry^ 

Taking the z-axis in the (001) direction (z will be con- 
sidered throughout this paper as the growth direction of 
the quantum well) and e = ez, the contribution of the 
Rashba term to the effective magnetic field may be writ- 
ten as 



Oii(k) = a ( -k 3 




(4) 



where a is the Rashba parameter. Next, we consider the 
influence on the effective magnetic field when we change 
the injection direction to (cos#, sin#, 0) and choose this 
direction as the cc-axis (A schematic of the configuration 
is shown in Fig. []}. The effective magnetic field changes 
according to fi'(k) = [/ _1 fi([/k), with 



cos 9 sin 9 ' 
V = I -sinfl cosfl 
1, 



(•5) 
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Therefore, we obtain for the Dresselhaus field: 
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(6) 



where the wave vector components are written in the ro- 
tated coordinate system. For the Rashba field we obtain: 



nk(k) = 




Notice that the Rashba field is expressed in the new co- 
ordinates in the same way as in Eq. (|rJ|). 

In the following, we discuss the diffusion in two differ- 
ent ways: (i) We take the identical strengths of a = (3 
with (3 = 7(^2) to study the diffusion direction depen- 
dence and the spin polarization dependence of the diffu- 
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sion. In this case, the spin orbit coupling £1' = Q' D + £l' R is written as 



n'(k) = 2/3 ($m(6 



cos 



n +7( — 



■ sin 29 + k x k y cos 29) 



(8) 



r 



(cos(7r/4 — 
sin(7r/4 — fj | repre- 


senting the crystal direction (110) in the rotated coordi- 
nates, (ii) We study the effect of the gate voltage on the 
spin diffusion. We treat this issue in Sec. IV C. 



The diagonal elements /k CT stand for the electron distri- 
bution functions of spin a whereas the off-diagonal ele- 
ments pktl = Pkit represent the correlations between the 
spin-up and -down states. 



III. KSBE AND INHOMOGENEOUS 
BROADENING 



Our fully microscopic treatment concerns the calcu- 
lation of the spin density matrix for the electron with 
momentum k at position r = (x,y): 



Pk(r, t) 



hi) 



(9) 



By using the non-equilibrium Green function method 
with gradient expansion as well as the generalized 
Kadanoff-Baym ansatz, 2 - we construct the KSBE as 
follows, 5 assuming the diffusion to take place in the x- 
direction: 



dp k (x,t) d^(x,t) dp k (x,t) 
h — h e- 



Ot 
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dk x 



h 2 k x dpk{x,t) 
to* dx 
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(10) 
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Here ^(x, t) is the electric potential satisfying the Poisson 
equation 



Vj*(r)=e[n(r)-JV (r)]/(oKo) 



(11) 



with n(r) standing for the electron density at position r, 
N (r) representing the background positive charge den- 
sity and kq being the static dielectric function. £ we (x,t) 
is the Hartree-Fock term from the Coulomb interac- 



tion. The scattering terms 



dt 



include the electron- 



electron and electron-phonon scattering. Their expres- 



sions can be found in Ref. |g. Details of the numerical 
scheme as well as the material parameters are laid out in 
Ref. [24]. It has been shown that Coulomb scattering plays 
a very important rule in spin dephasing & 21 i 25 i 26 and spin 
diffusion/transport— and therefore cannot be neglected 
in calculating the spin diffusion length. 



When the system reaches its steady state, in the ab- 
sence of an external field, it can be described by the equa- 
tion 



J 



h 2 k x dp k (x,t) 
m* dx 
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(12) 



By neglecting the Hartree-Fock and the scattering terms, Eq. (TT2")) can be solved analytically. 24 After dividing both 
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sides of the equation by k x , we see that along with the 
diffusion, the spin polarization for each wave vector k pro- 
cesses along wit = 2 ^ k Sl'(k). The fact that spins with 
different momentums precess with different frequencies is 



referred to as inhomogeneous broadening] i 28 With any 
spin-conserving scattering included, the inhomogeneous 
broadening results in spin dephasing.^21 
For a — (3, reads 



m 
2ft 2 
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(13) 



We find that u>k includes two parts: the zeroth-order 
term (on k) which is always along the same direc- 
tion fio, and the second-order term. The zeroth- 
order term includes both the k-dependent and the k- 
independent terms, while the second-order term is always 
k-dependent. According to the previous workS ) 5 ' 6 ' 24 i 27 
the k-dependent term contributes to the inhomogeneous 
broadening. 

We hrst analyze the spin diffusion without the third- 
order term of the Dresselhaus term (and hence the 
second-order term of ut). In this case, the effective mag- 
netic field for each wave vector k points to the same di- 
rection no. From the previous works ) 10 i 21 we know that 
the spin dephasing in the time domain, in the spatial ho- 
mogeneous case, shows strong anisotropy with respect to 
the direction of spin polarization: The spin polarization 
along the direction no has an infinite dephasing time, 
while the ones vertical to this direction have a very short 
spin dephasing time due to the large spin-orbit coupling. 
This is because when the spin polarization direction is 
the same as the direction of the effective magnetic field, 
the spin polarization cannot precess and hence there is 
no inhomogeneous broadening. Consequently the spin 
polarization has an infinite dephasing time. The same is 
true for the spin diffusion: When the spin polarization 
is along the same direction of w^, i.e., fio (the crystal 
direction (110)), the spin polarization will not decay re- 
gardless of the direction of the spin injection, even in the 
presence of scattering, as there is no spin precession and 
hence no inhomogeneous broadening; However, when the 
spin polarization is perpendicular to no, the spin diffu- 
sion length is very short due to the large inhomogeneous 
broadening. 

Nevertheless, in contrast to the spin dephasing in a 
spatially homogeneous system, the spin diffusion has an 
extra degree of freedom, i.e., the spin diffusion/injection 
direction (which should not be confused with the spin 
polarization direction). This degree of freedom intro- 
duces an additional level of anisotropy to the spin dif- 
fusion/transport. When the spin diffusion/injection di- 
rection is along 8 = 3tt/4, i.e., (110) direction, the pre- 
cession frequency = ^^2^n becomes k-independent 
if the cubic term of the Dresselhaus term is neglected. 
Moreover, even with the scattering included, Eq. (TT2")) 



can also be satisfied by the following solution: 

p w ( X ) = e ^-y{k 2 z )a v x/h' < p ^ x = ty e -irn'-y{k 2 z )* v x/h 2 ^ 

where pk(x = 0) is the Fermi distribution with the spin 
polarization at the left boundary. We see from Eq. (fl4|) 
that in this special case the oscillation period is the same 
for all k's, with the value — , ,, ■>, . So, there is no inhomo- 
geneous broadening and the spin has an infinite diffusion 
length, regardless of the direction of the initial spin po- 
larization. This result is mainly induced by the peculiar 
precession frequency which is determined by the spin- 
orbit coupling and the diffusion velocity k x , so it cannot 
be obtained in the spatially homogeneous case where the 
precession frequency is only determined by the spin-orbit 
coupling. Of course, when we include the cubic term in 
the spin-orbit coupling, the oscillation period becomes 
wave- vector dependent. This inhomogeneous broaden- 
ing, then, results in a finite diffusion length. 

In brief, we stress the importance of the inhomoge- 
neous broadening caused by the diffusion term as com- 
pared to the spin dephasing in the spatially homogeneous 
case. Consequently the spin polarization shows strong 
anisotropy not only for the direction of the initial spin 
polarization but also for the direction of the spin diffu- 
sion. Moreover, the cubic Dresselhaus term cannot be 
neglected. In what follows we solve the problem by using 
the fully microscopic many-body approach of the KSBE 
with the cubic Dresselhaus term and the scattering ex- 
plicitly included. 



IV. NUMERICAL RESULTS 

The KSBE are solved numerically for (001) GaAs 
quantum wells of width a = 5 nm including the electron- 
electron and the electron-LO-phonon scattering. The 
value of the Dresselhaus coefficient is as 7 = 25 eV-A 3 ^ 
The electron density and the temperature are taken as 
N e = 4.0 x 10 11 cm" 2 and T — 200 K separately unless 
otherwise specified. All matrix elements of the interac- 
tions are given in Ref. |8|. For the spin polarization along 
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FIG. 2: Spin diffusion length Ld (solid curves) and the inverse 
of the spin oscillation period Lq (dashed curves) for a = /3 
as functions of the injection direction for different spin polar- 
ization directions no, z and ni at T = 200 K. It is noted that 
the scale of the spin oscillation period is on the right hand 
side of the frame. 



direction n, the boundary conditions are given by 21 




o>*)L>o = 

L,t)\ k 



^k.t — Fir 



-n • a 



k x <0 



(15) 

with the Fermi distribution F Ka = [ e (fc 2 /2m*- M „)/fc B T + 
{a =tl) and the chemical potential \x a determined 
by the polarized electron density. The diffusion length 
and the oscillation period are extracted from the spa- 
tial evolution of the spin polarization along n direction 
AN(x) = J2k 1' r [Pk(a;)n • <x]. 



A. Spin-diffusion/injection— direction and 
spin-polarization— direction dependence at a = 

We first fix the spin polarization along the (110) direc- 
tion (no) and study the spin diffusion length as a function 
of the spin injection direction in the presence of the cu- 
bic Dresselhaus term. For comparison, we also study the 
cases when the spin polarizations are perpendicular to 
the (110) direction, i.e., z and ni = z x no and show 
how the spin injection lengths change as a function of 
spin diffusion/injection direction. Subsequently we ex- 
plore the special case when the spin diffusion/injection 
direction is along (110), i.e., 9 — 37r/4 and show how the 
spin injection length changes with the spin polarization 
in the presence of the cubic Dresselhaus term. The spin 
polarization calculated from the KSBE Eq. (JTUJ) and the 
Poisson Eq. (fTTj) along these three spin polarizations can 
be well fitted by 



AN(x) = Cexp(-^-)cos(^ 



(16) 



where Ld is the spin diffusion length and Lq represents 
the oscillation period. 

In Fig. [5] the spin diffusion length Ld and the inverse 
of the spin oscillation period Lq 1 are plotted against 
the spin diffusion/injection angle 9 for spin polarizations 
along n , z and ni, respectively. It is interesting to see 
that with the inclusion of the cubic Dresselhaus term, 
the spin injection length becomes finite but still inde- 
pendent on the direction of the spin diffusion/injection if 
the spin polarization is along (110) (= no). This is in 
agreement with the case without the cubic term where 
Ld becomes infinite. At the same time, the spin po- 
larization along this direction has an infinite oscillation 
period Lq = oo. This can be understood because the 
spin polarization is in the same direction as the effective 
magnetic field given by the zeroth-order term in u>k and, 
thus, results in no oscillations. It is noted that the ef- 
fective magnetic field from the second-order term in Wk 
is k-dependent and therefore cannot lead to any oscilla- 
tion at high temperature due to the scattering . 21 ' 25 ' 30 
However, this second-order term causes an inhomoge- 
neous broadening which leads to a finite spin diffusion 
length. By rewriting the second-order term of into 

, 2 (ky/kx\ 
2p^sin(20 k + 26>) -1 with k and (9 k denoting 

v /. 

the magnitude and the direction of the wave vector k 
separately, one finds that the magnitude of the inho- 
mogeneous broadening does not change with the spin 
diffusion/injection direction 9. Consequently the spin 
diffusion length does not change with the spin diffu- 
sion/injection direction. 

Strong anisotropy is again found for spin polarizations 
along the directions perpendicular to n except for the 
special case when 9 = 3ir/4. Moreover, it is found in 
the figure that, for the perpendicular directions, both Ld 
and Lq depend sensitively on the spin diffusion/injection 
direction (additional anisotropy). It is further noticed 
that Ld and Lq are almost identical for both the per- 
pendicular directions z and h\. These results can be 
understood again from the inhomogeneous broadening 
induced by u>k- As 7 in Eq. (|13j) is much smaller than 
(3, the inhomogeneous broadening for spin polarizations 
z and ni are therefore defined by the zero-th order term 
of u>k, which is the same for both perpendicular direc- 
tions and is spin diffusion/injection direction 9 sensi- 
tive. The oscillation period in the diffusion Lq is de- 
termined by the k-independent component of Wkr^ i.e., 
m* /3sin(6> — ir/A)/h 2 which is in good agreement with the 
results in Fig. [2] 

Now we turn to the case with special fixed spin dif- 
fusion/injection direction 9 = 3tt /4, and investigate the 
spin diffusion length with different spin polarizations. It 
is interesting to see from Fig. [5] that similar to the case 
without the cubic Dresselhaus term, the spin diffusion 
lengths for the three spin polarizations are almost iden- 
tical. Therefore the anisotropy disappears for this par- 
ticular spin diffusion/injection direction. This is because 
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FIG. 3: Spin diffusion length Ld and spin oscillation period 
Lo at 9 = 37r/4 and T — 200 K as functions of the electron 
density. Note the scale of the oscillation period is on the right 
hand side of the frame. 
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FIG. 4: Spin diffusive length Ld and spin oscillation period 
io at two different electron densities, 4 x 10 11 cm~ 2 and 4 x 
10 10 cm -2 vs. the temperature. The spin diffusion/injection 
direction 8 = 3n/4. Note the scale of the oscillation period is 
on the right hand side of the frame. 



sion length are plotted as functions of the electron den- 
sity. One finds that the diffusion length Ld decreases 
with the density, whereas the oscillation period Lo in- 
creases with it. However, it is noted that the change 
observed in the oscillation period is almost negligible (1 
%) in contrast to the pronounced changes observed in the 
diffusion length (of the order of 40 %). We further inves- 
tigate the temperature dependence of the spin diffusion 
length and the oscillation period in Fig. [4] by taking two 
typical electron densities, i.e., N — 4 x 10 10 cm~ 2 and 
4 x 10 11 cm~ 2 . Similar behavior is found here compared 
to the density dependence: the spin diffusion length Ld 
decreases markedly with the temperature whereas the os- 
cillation period Lo only increases slightly with it. 

Both dependences above can be again understood by 
the inhomogeneous broadening. Unlike the results in Rcf. 
[24j where only the Dresselhaus term is considered and a 
very weak temperature dependence is obtained, the dif- 
fusion length here has a strong temperature dependence. 
This can be understood by the difference in the behavior 
of the inhomogeneous broadening, which is determined 
by the anisotropy of the precession frequency ut. In the 
previous case, the zero-th order term dominates the in- 
homogeneous broadening which changes little with tem- 
perature and electron density. However, in the present 
case with a = f3 and 9 = 3n/4, the k-dependent zeroth- 
order term is always zero and the second-order term alone 
determines the inhomogeneous broadening. This term 
increases effectively with the temperature and with the 
electron density due to the increase of the average value 
of k 2 . This is the reason for the obtained marked de- 
crease of the spin diffusion length. The oscillation pe- 
riod Lo is determined by the k-independent zeroth-order 
term of ut, which is independent of the electron den- 
sity and the temperature. Therefore one observes only a 
slight change in the oscillation period, originating from 
the second-order term and the scattering. 



Gate-voltage dependence with injection 
direction along (110) 



when the spin diffusion/injection direction is 9 = Sir/4,, 
the k-dependent component in the zero-th order term of 
u>k disappears. This results in the same inhomogeneous 
broadening for any spin polarization. 



B. Temperature and electron-density dependence 
with injection direction along (110) and a — f3 

Now we investigate the temperature and electron den- 
sity dependence of the spin diffusion length at the special 
spin diffusion/injection direction 9 — Sir/A with a = (3. 
As shown in the previous section the spin diffusion is 
isotropic with respect to the spin polarization direction. 
Hence we choose the spin polarization to be along the 
z-axis. In Fig. [3] the oscillation period and the diffu- 



We now tune the gate voltage in order to change the 
relative importance of the Dresselhaus and the Rashba 
terms. The spin diffusion/injection direction is again 
fixed at 9 — 37r/4. The calculated spin diffusion length 
Ld and the spin oscillation length Lo are plotted against 
a//3 in Fig. O at T = 200 K, with the spin polarizations 
being along no, z and ni respectively. It is noted that 
since Lo = oo when the spin polarization is along no 
[(110)], we only plot Lo for the two perpendicular direc- 
tions in the figure. 

Three main results are obtained: (i) Differing from the 
case without the cubic Dresselhaus term, the maximum 
diffusion length occurs interestingly at a/f3 ~ 0.95 (in- 
stead of 1) for each spin polarization direction; (ii) The 
isotropy of the spin polarization direction at a//3 =/= 1 is 
suppressed except at a/ (3 ~ 0.9, and the diffusion lengths 
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FIG. 5: Spin diffusion length Ld (solid curves) and spin os- 
cillation period Lo (dashed curves) as functions of the ratio 
a/f3 with the injection direction 9 = 3tt/4 for different spin 
polarization directions no, z and fii at T = 200 K. The chain 
line is used to guide the eyes. The scale of the oscillation 
period is on the right hand side of the frame. 



for the spin polarization along z and fii are always iden- 
tical; (iii) The spin oscillation period Lq decreases with 
a/ (3. These results can also be understood from the in- 
homogeneous broadening. For 9 — 37r/4, the spin pre- 
cession frequency along the diffusion direction is given 
by 



[R3 + a+^)^-7^] ni 
- [(3 + a - l[kl - k 2 y )]h . 



(17) 



At T = 200 K, (fc 2 )/2 ~ 0.11(7r/a) 2 with (•) representing 
the average over the imbalance of the spin-up and -down 
electrons. Therefore, (—/3 + a + ^^-) ^-fii is almost zero 
at a/ (3 ~ 0.95 and does not contribute to the inhomo- 
geneous broadening. This results in a maximum of the 
spin diffusion length. It is further noted that the inho- 
mogeneous broadening term ^(/c 2 — fc 2 )no gives a larger 
effect when the spin polarization is along z and ni than 
that when the spin polarization is along n as the lat- 
ter coincides with the precession axis. This explains why 
one obtains a larger diffusion length when the spin po- 
larization is along no at a/ (3 ~ 0.95. Finally as the spin 
oscillation period L is determined by the magnitude of 
the k-independent term of Uk, i-e., (3(a/(3 + l)n , which 
increases with a/ (3. Therefore, L decreases with a//?. 



V. CONCLUSIONS 

In conclusion, wc study the spin diffusion in (001) 
GaAs quantum wells with competing Dresselhaus and 



Rashba spin-orbit coupling strengthes by solving the 
KSBE with the electron-phonon and the electron- 
electron Coulomb scattering explicitly included. It is 
shown that unlike the spin dephasing in the time do- 
main where strong anisotropy is determined by different 
spin polarization, for spin diffusion it is also determined 
by the spin diffusion/injection direction. By neglecting 
the cubic term of the Dresselhaus spin-orbit coupling and 
with a = (3, the ideal case of an infinite diffusion length 
is obtained either for the spin polarization along (110) 
regardless of the spin diffusion direction or for the spin 
diffusion direction along (110) regardless of the spin po- 
larization direction. However, the cubic term cannot be 
neglected, resulting in a finite diffusion length, which in 
fact is small, about 2 /im at the temperature of 200 K 
when the electron density lies between 4 x 10 10 cm~ 2 and 
4 x 10 11 cm~ 2 . It is then shown that when a = (3, for the 
spin polarization along (110), the spin diffusion length 
changes little for the spin diffusion direction, whereas for 
the other two perpendicular spin polarization directions 
(110) and (001), the spin diffusion length shows strong 
anisotropy for the spin diffusion/injection direction and 
has a peak when the spin diffusion/injection direction is 
along (110). When the spin diffusion/injection direction 
is along (110), the spin diffusion length is isotropic for 
spin polarization. The electron density and temperature 
dependence of the spin diffusion is closely investigated. 

We also tune the gate voltage to show the competing 
effect of the Dresselhaus and the Rashba strengthes. It is 
found that with the cubic Dresselhaus term included, the 
maximum spin diffusion length appears at a/ (3 ~ 0.95 in- 
stead of 1 where the cubic term is neglected. The spin dif- 
fusion length decreases when the ratio of the Rashba and 
Dresselhaus strengths deviates from 0.95. These results 
can be well understood by the inhomogeneous broaden- 
ing. 
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